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Abstract 

We report the observation of a new D s j meson produced in B + — > D°D s j — > D°D°K + . This 
state has a mass of M = 2715 ± llt^ MeV/c 2 , a width r = 115 ± 20+32 MeV/c 2 and a spin-parity 
1~. The results are based on an analysis of 449 million BB events collected at the T(4S I ) resonance 
in the Belle detector at the KEKB asymmetric energy e + e~ collider. 

PACS numbers: 
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At the level of quark diagrams, the decay B — > DDK proceeds dominantly via the 
CKM-favored b — > cW + — > ccs transition. The transition amplitudes can be categorized 
as either external W - or internal W -emission (color-suppressed) diagrams. The decay 
B + — > D°D°K + 1] can proceed through both types of diagrams, with the naive expectation 
that the internal W contribution to the branching fraction is suppressed relative to that 
of the external W by a factor of nine. 5-meson decays to three-body DDK final states 
are a promising area for searches for new cs states as well as of some cc states lying above 
D°D° threshold. Since the externally emitted W produces 1 + , 1~, CT states, cs mesons with 
these quantum numbers should be copiously produced. The unexpected discoveries of the 
D s j(2317) and -D s j(2457) mesons show that our understanding of the cs spectroscopy might 
be incomplete. Experimental data on cc states with decay channels open to D^D^ are 
scarce. 

The decays B — > DDK have been previously studied with a small data sample by 
ALEPH 0. Recently a comprehensive study was performed by BaBar 0. In this letter we 
report the first study of the Dalitz plot of B + — ► D°D°K + decay. 

The study is performed using data collected with the Belle detector at the KEKB 
asymmetric-energy e + e~ (3.5 on 8 GeV) collider j^|, operating at the T(4S) resonance 
(-y/s = 10.58 GeV) with a peak luminosity that exceeds 1.6 x 10 34 cm~ 2 s _1 . The data 
sample corresponds to the integrated luminosity of 414 fb _1 and contains 449 million BB 
pairs. 

The Belle detector is a large- solid- angle magnetic spectrometer that consists of a sili- 
con vertex detector, a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons. 
The detector is described in detail elsewhere (if. 

We select charged tracks that originate from the interaction region by requiring | dr | < 
0.4 cm and |dz|< 5 cm, where dr and dz are the distances of closest approach to the inter- 
action point in the plane perpendicular to the beam and along the beam axis, respectively. 
Charged particles are identified by using combined information from the TOF, ACC and 
dE/dx measurements in the CDC. Requirements on the particle identification variable are 
imposed that identify a charge kaon with 90% efficiency at < 10% tt —>■ K miss-identification 
probability and a charged pion with « 100% efficiency at < 10% K — > tt miss-identification 
probability. Any track that is positively identified as an electron is rejected. 

Candidate Kg — > tt + tt~ decays are identified by a displaced secondary vertex, a two-pion 
momentum vector that is consistent with a Kg originating from the interaction point and 
by the invariant mass selection | M n + n - — m K o |< 15 MeV/c 2 . Candidate tt° mesons are 
identified as pairs of ECL-identified photons, each with a minimum energy of 50 MeV, that 
have an invariant mass within ±15 MeV/c 2 of the 7T° mass. 

D° mesons are reconstructed in the K~tt + , K~tt + it + it~ , K~tt + it , i^°7r + 7r~ and K~~K + 
decay modes. We preselect D candidates using a signal window ±30 MeV/c 2 around the 
nominal D meson mass for decay modes, except for D° —>■ K~tt + it decays, where a larger 
±50 MeV/c 2 signal window is used. Mass- and vertex-constrained fits are applied to all D 
meson candidates to improve their momentum resolution. 

We retain events that have a kaon candidate and at least two candidate D° mesons, 
with allowed flavor combination and with the D° momenta in the T(4S) rest frame ('cms') 
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below the kinematical limit for B + — > D°D°K + . To suppress the continuum background 
(e + e~ — > qg, q — u,d, s, c) we require the ratio of the second to the zeroth Fox- Wolfram 
moments 6;] to be less than 0.3. Momenta of the secondaries from a B meson candidate 
decay are refitted to a common vertex with an interaction point constraint that takes into 
account the B meson decay length. 

The B meson candidates are identified by their cms energy difference, AE = Eji?j— i?beam> 
and their beam constrained mass, M bc = \J E^ eam — (Sj^) 2 , where E heSLm = y/s/2 is the beam 
energy in the cms and Pi and Ei are the three-momenta and energies of the B candidate's 
decay products. For the subsequent analysis we select B candidates with M bc > 5.2 GeV/c 2 
and -0.4 GeV < AE < 0.3 GeV. Exclusively reconstructed B + — > D°D°K + signal events 
have Mb c distributions that peak at the nominal I?-meson mass; the AE distributions peak 
near zero. 

We employ a discriminator (likelihood ratio) based on the D° meson signal significance 
to select the unique B candidate in the event, defined as: C7Z(Md) = s ^[ S ^+b\m d ) > wnere 
S and B are the signal and the background probabilities respectively. This discriminator 
is determined from the data, for each D° decay mode separately, using a sample enriched 
in B — > DDK decays. For events with multiple B — > DDK candidates, the product 
CIZb = CTZ(Mjj) x C1Z(Md) is calculated and the candidate with the largest CTZ(B) is 
accepted. The solution with larger kaon identification likelihood C(K) is chosen if multiple 
kaon candidates are found accompanying the accepted D°D° combination. 

The CIZb discriminator is also used to suppress combinatoric backgrounds to B + — > 
D°D°K + and to enhance the signal purity. 

The AE and Mb c distributions for the B + — > D°D°K + decay candidates, selected with 
CIZb > 0.01 requirement, are shown in Fig. ^ The AE distribution is shown for events 




with | M bc — m B \< 3(T A ,/ bc (cr Mhc = 2.7 MeV/c 2 ), where m B is the nominal B meson mass 
and the M bc distribution is shown for events with | AE \ < 3<ja_e (o"ab = 6.6 MeV). 

From a study of the M^ c and AE background distributions in large MC samples of generic 
BB and qq events, we find that the background level in the data is well explained by the MC 
simulations and no peaking component is seen in either distribution. Using events where 
one or both D° candidates are from the D°-mass sidebands, we verify with the data that 
there is no significnat peaking background. 

To extract the signal yield, we perform extended unbinned maximum-likelihood fits si- 
multaneously to AE and Mb c . The probability density functions (PDFs) for the Mb c and 



5 



AE signals are Gaussians. The background PDF for M^ c is represented by a phenomeno- 
logical function 



with a phase-space-like behaviour near the kinematic boundary; the AE 
background is parameterized by a second-order polynomial. The likelihood function is max- 
imized with free parameters for the signal yield, the Gaussian means and widths, and four 
parameters that describe shapes of the background distributions. 

From the fit, we obtain a signal yield of iV si g = 399 ±40 events with a signal-to-background 
ratio of S/B ~ 0.3. The results of the fit are superimposed on the AE and M^ c projections 
shown in Fig. ^ We determine the branching fraction from the relation 
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where are efficiences for the i and j D subchannels and for N B + B - ; we assume N B + B - = 
N B ogo. The efficiences are determined by MC using a phase-space decay model. The sum 
in the denominator of Eq.Q amounts to 6.8 x 10 -4 . 

We obtain B(B + -> D°D°K + ) = (13.1 ±1.3±$) x 10" 4 , where the first error is statistical 
and the second is systematic. The latter includes contributions due to uncertainties in the 
efficiency determination (tracking and particle identification efficiency, data-MC differences 
in AE, Mb c signal shapes), the CTZ B selection, the background parameterization, the MC 
model used in the efficiency calculation, the intermediate D — > i branching fractions and 
N B+B -. 
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FIG. 2: Dalitz plot (a) and projections for 1? + — ► D°D°K + in 1.5<7 AE — M^ c signal region: 
M(D°D°)(h) : M(D°K + )(c), M(D°K + )(d). Hatched histograms represent background, red curves 
show the efficiency variation (vertical scale in %). 

The Dalitz plot M 2 (D°D°) vs M 2 (D°K + ) for events from a signal region defined by 
the ellipse (AE/1.5a AE ) 2 + ((M bc - m B )/1.5a Mbc ) 2 < 1 is shown in Fig.^a). The three, 
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two-body invariant mass distributions are shown in Figs. |2(b)-(d). The hatched histograms 
represent the background distributions obtained for events from an elliptical strip surround- 
ing the AE, Mb c signal region that extends from 6<7AE,0M bc to 10a ae, &M bc - The back- 
ground distributions are normalized to the number of background events under the signal 
peaks (±1.5<r) as determined from the combined AE and M^ c fit. The data shown are not 
efficiency corrected. The variation of efficiency as a function of invariant mass is shown in 
Figs. 12(b)- (d) as a continuous curve. 

A pronounced feature of the Dalitz plot is the accumulation of events in the region 
16 GeV 2 /c 4 < M 2 (D°D°) < 18 GeV 2 /c 4 and 7 GeV 2 /c 4 < M 2 (D°K+) < 8 GeV 2 /c 4 , possibly 
the overlap of a horizontal band that could be due to the 0(4160) and a vertical band 
that cannot be attributed to any known cs state. A horizontal band at M 2 (D°D°) ~ 
14.2 GeV 2 /c 4 corresponds to the ■0(3770) production. 

We employ simultaneous fits to the AE and M^ c distributions for events from each 
50 MeV/c 2 mass bin of the Dalitz plot projection to obtain background-subtracted invariant 
mass distributions. In these fits the widths and positions of the Gaussians describing the 
signal are fixed at the values obtained for the total signal sample, while the signal yield 
and the background PDF's paremeters are free parameters. The obtained signal yields are 
shown in Fig. 0] as points with error bars. 

The '0(3770) signal is studied with finer, 20 MeV/c 2 mass bins (Fig. 0(a)). The peak 
is fitted in the region M(D°D°) < 4 GeV/c 2 with a Breit-Wigner (BW) plus a threshold 
function to describe a nonresonant component. The 0(3770) signal yield is 68 ± 15 events 
with a peak mass of 3777 ± 3 MeV/c 2 , and a width of 27 ± 9 MeV/c 2 , in agreement with 
the PDG averages. 
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FIG. 3: Background-subtracted mass distributions: (a) -0(3770) signal, (b) M(D°D°) for cosOhel > 
(c) M(D°K + ) for M(D°D°) > 3.85 GeV/c 2 Solid curves denote fit results described in the text. 
The lower curve in (a) shows phase-space component, coloured curves in (c): purple - the tp(4:l60) 
reflection, green: phase-space, red: threshold component and navy-blue: the sum of the three 
components. 

The background-subtracted M(D°D°) spectrum (Fig. E(b)), for events satisfying 
cosOhel > 0, where 6h e i is the helicity angle between the D° momentum vector and the boost 
direction to the D°D° rest frame, is used to estimate the possible 0(4160) contribution to 
the enhancement at M(D°K + ) ~ 2.7 GeV/c 2 . The distribution for M(D°D°) > 3.8 GeV/c 2 
is fitted with a BW with mass and width fixed at the nominal ^(4160) values (M = 4160, 
T = 80 MeV/c 2 j8j), yielding 24 ± 11 signal events. We use these 0(4160) parameters to 
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estimate the number of -0(4160) events in the backward helicity-angle hemisphere, in the 
region M{D°K + ) < 2.9 GeV/c 2 . 

Figure Efc) shows the background-subtracted M(D°K + ) distribution for events with 
M(D°D°) > 3.85 GeV/c 2 . This requirement removes the -0(3770) reflection at high 
M(D°K + ). The predicted -0(4160) reflection is indicated in the figure by the purple curve. 
The -0(4160) reflection agrees well with the data in the high mass M(D°K + ) region but 
does not explain the large peak at M(D°K + ) ~ 2.7 GeV/c 2 . We parameterize the observed 
excess of events with a BW and fit the M(D°K + ) spectrum (Fig. Etc)) with the ansatz of 
a new resonance, the 0^(4160) reflection shape and a phase-space component as determined 
by MC simulations. The free parameters in the fit are the resonance yield, mass and width, 
and the phase-space component normalization. The fit has an acceptable overall x 2 but 
is unable to reproduce the events near the low- mass threshold seen in Fig. E^c). We used 
several phenomenological parameterizations (polynomials, another BW, an exponential) of 
the threshold enhancement into the fit to determine its influence on the BW parameters of 
the 2.7 GeV/c 2 peak. The exponential form a x exp[— aM 2 (D°K + )] gives a good description 
of the mass spectrum, while adding only two free parameters. From this fit we obtain for 
this new resonance, which we further denote as the _Dj}(2700), a signal yield of 182 ± 30 
events, mass of M = 2715 ±11 MeV/c 2 and width of T = 115 ± 20 MeV/c 2 . The threshold 
and the phase-space components from the fits are 58 ± 38 and 47 ± 26 events, respectively. 
The fit results are depicted in Figs. E^a)-(c) as histograms overlaid on the measured mass 
spectra. 

The resonance parameters and product branching fractions are summarized in Table 
The -0(4160) yield is not statistically significant, and therefore a 90% C.L. upper limit is also 
quoted. The £>'s of the threshold and the phase-space components are (1.9±1.2±£;S) x 10~ 4 
and (1.5 ± 0.8i?;g) x 10 -4 (the first errors are statistical, the second are systematic) which 
correspond to the 90% C.L. upper limits of 4.9 x 10~ 4 and 3.0 x 10~ 4 , respectively. The 
systematic errors on the product branching fractions and the resonance parameters include 
contributions from the efficiency variation over the Dalitz plot, uncertaintities in the yields 
of the -0(4160) reflection (including recent measurements of the -0(4160) parameters [9]), the 
threshold parameterization, sensitivities of parameters to the fit range and parameteriza- 
tion, uncertainties in the CIZb selection, as well as due to neglected interference effects. The 
systematics due to the latter are determined from MC simulations of Dalitz plot densities 
with and without interference of coherent amplitudes. Here an isobar formalism was used, 
with each contributing resonance parameterized by the BW form. The resonance parameters 
from Table |l] and the threshold enhancement parameters are used to determine the ampli- 
tudes. The effects of interference of the -0(3770) with other states are found small and are 
neglected in the simulations. Events generated with maximal constructive and destructive 
interferences between the amplitudes, were passed through a complete detector simulation 
and analysed ignoring interference effects. The differences between the obtained resonance 
parameters and the input values are taken as their systematic errors. 

We study background-subtracted 0(3770) and Z)jj(2700) helicity angle distributions by 
selecting the respective invariant mass in the resonance region and obtaining signal yields 
in bins of cos 6^/ from simultaneous fits to AE and Mt, c . Spin-parity hypotheses for the 
resonances are tested by performing binned x 2 fits to the obtained angular distributions 
corrected with efficiency weights. The J = 1 hypothesis describes the 0(3770) distribution 
well (x 2 /ndf = 3.6/5), J = 1 is favoured for the D+j(2700) (jf/ndf = 7/5); the J = 
(185/5) and J = 2 (250/5) assignments can be rejected. 



TABLE I: Branching fractions of quasi-two-body components 
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FIG. 4: 5+ D°D°K + signal yield vs M(D°D°) (a), M( J D°K+) (b), and M(D°K + ) (c). 
Histograms denote the contributions from: Dj^(2700) (blue), "0(3770) (green), 0(4160) (yellow), 
threshold (red) and phase-space components (navy blue). Histograms are superimposed additively. 



The M(D°K+) distribution of Fig.EJc) in 10 MeV/c 2 is used to search for the L> sJ (2573) 
contribution. We include this state in the fit using the BW formula with the parameters 
from j.8], and obtain N sig (D s j (2573)) = 7.7 ± 5.7, which corresponds to the 90% C.L. upper 
limit of B{B+ -> D°D sJ {2573)) x B{D sJ (2573) -> D°K+) < 0.7 x 10~ 4 . 

In summary, from a study of the Dahtz plot we find that the decay B + — > D°D°K + pro- 
ceeds dominantly via quasi-two-body channels: B + — > D°Dfj(2700) and B + — > ip(3770)K + , 
where Djj(2700) is a previously unobserved resonance in the D°K + system with a mass 
M = 2715 ± lltll MeV/c 2 , width T = 115 ± 20ijj!j MeV/c 2 and J p = 1". The observed 
rate for -0(3770) production confirms our previous observation 

Based on its observed decay channel, we interpret the Dfj(2700) resonance as a cs meson. 
Potential model calculations predict a cs radially excited 2 3 5*i state with mass M = 
2720 MeV/c 2 . From chiral symmetry considerations [l2| a 1 + -1 _ doublet of states has been 
predicted. If the 1 + state is identified as the _D s j(2536), the mass predicted for the 1~ state 
is M = 2721 ± 10 MeV/c 2 . Additional measurements of the meson properties are needed to 
distinguish between these two interpretations. 
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